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Abstract Inspired by the squared eigenfunction symmetry constraint, we intro- 
duce a new T^-flow by "extending" a specific £„-flow of discrete KP hierarchy 
(DKPH). We construct extended discrete KPH (exDKPH), which consists of t n - 
flow, Tfc-flow and t n evolution of eigenfunction and adjoint eigenfunctions, and 
its Lax representation. The exDKPH contains two types of discrete KP equation 
with self-consistent sources (DKPESCS). Two reductions of exDKPH are ob- 
. tained. The generalized dressing approach for solving the exDKPH is proposed 

\ and the N-soliton solutions of two types of the DKPESCS are presented. 

• i-H ■ 

1 Introduction 

Generalizations of soliton hierarchy attract a lot of interests from both phys- 
ical and mathematical points and there were some methods to generalize the 
soliton hierarchy[TJ-[3]. Recently, a systematic approach inspired by squared 
, eigenfunction symmetry constraint was proposed to construct the extended KP 

£NJ ■ hierarchy [5]. By this method, the extended two-dimensional Toda lattice hier- 

j — [ , archy, the extended CKP hierarchy and the extended q-deformed KP hierarchy 

<^ ' have been obtained[6]-[8]. 

ON . The discrete KP hierarchy (DKPH) [9]- [12] is an interesting object in the 

research of the discrete integrable systems and the discretization of the inte- 
grable systems [15]. The Sato's approach for the discrete KPH was presented 
in [TT] . Naturally, there are some similar properties between discrete KPH and 
KPH[]3], such as tau function [T2l [T^] . Hamiltonian structure[H] and gauge 
transformation[10l [T51 [16], etc. In [10], Oevel has given explicitly two types 
of gauge transformation operators of the discrete KPH. In [TB], the combined 
gauge operator and the determinant representation of the operator have been 
obtained. 

In this paper, we will construct the extension of the discrete KPH(exDKPH). 
Inspired by the squared eigenfunction symmetry constraint of discrete KP hi- 
erarchy [10], we introduced the new Tfc-flow by "extending" a specific i„-flow 
of discrete KP hierarchy. Then we find the exDKPH consisting of i„-flow of 
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discrete KP hierarchy, Tfc-flow and the Revolutions of eigenfunctions and ad- 
joint eigenfunctions. The commutativity of i n -flow and r^-flow gives rise to zero 
curvature representation for exDKPH. Also the Lax representation of cxDKPH 
is derived. Due to the introduction of r^-flow the exDKPH contains two time 
series {t n } and {t^} and more components by adding eigenfunctions and ad- 
joint eigenfunctions. The exDKPH contains the first type and second type of 
discrete KP equation with self-consistent sources(DKPESCS). The KP equation 
with self-consistent sources arose in some physical models describing the inter- 
action of long and short waves [3]. The similarity of KP equation and discrete 
KP equation enables us to speculate on the potential application of discrete 
KP equation with self-consistent sources. By i„-reduction and r^-reduction, 
the exDKPH reduces to a discrete 1+1-dimcnsional integrable hierarchy with 
self-consistent sources and constrained discrete KP hierarchy, respectively. 

The dressing method is an important tool for solving soliton hierarchy [12] . 
However this method can not be applied directly for solving the "extended" 
hierarchy A generalized dressing approach for exKPH is proposed in |17j . In 
this paper, with the combination of dressing method and variation of constants 
method, a generalization to the dressing method for exDKPH is presented, which 
is based on the dressing method for discrete KPH [TT] and the similar approach 
for finding Wronskian solutions to constrained KP hierarchy [TH]. In this way, 
we can solve the entire hierarchy of exDKPH in an unified and simple manner. 
As the special cases, the N-soliton solutions of the both types of DKPESCS are 
obtained simultaneously. 

This paper will be organized as follows. In Sec. 2, we present the exDKPH 
and its Lax pair, which includes two types of DKPESCS. In Sec. 3, i„-reduction 
and Tfc- reduction for the exDKPH are given. In Sec. 4, we discuss the generalized 
dressing method for the exDKPH. In Sec. 5, we present the N-soliton solutions 
of the DKPESCS. 

2 New extended discrete KP hierarchy 

We denote the shift and the difference operators acting on the associative ring 
F of functions by T and A, respectively, as follows 

F = {f(l) = f(l,t 1 ,t 2 ,--- ,t i ,---);l€Z, t % eR} 

r(/(i)) = /(/ + 1) = / (1) (0, A(/(I)) = f(l + 1) - f(i). 

In this paper, we use P(f) to denote an action of difference operator P on the 
function /, while Pf means the multiplication of difference operator P and zero 
order difference operator /. Define the following operation 

Also, wc define the adjoint operator to the A operator by A* 

a*(/(0) = (r- 1 - J)(/(0) = f(i i) - f(i), (2) 
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A *'f = E ( l) ( A "^ + 1 - 3))) A * J ~ l - (3) 

Let P = Ej=-oo /j(0 Ai , thc adjoint operator P* is denned by P* = £j=-oo &* j fj(l)- 
The Lax equation of the DKP hierarchy is given bv[9l [TT] 

L tn = [B„,i], (4) 

where L = A + ,/o + ./iA~ 1 + /2A~ 2 + --- is a pseudo-difference operator with 
potential functions fi <E F, B n = I7\ stands for the difference part of L n . Thc 
commutativity of t n - and i m -flow gives rise to the zero-curvature equations for 
DKP hierarchy: 

B n ,t m — B m j n + [B ni B m ] = 0. (5) 
with the Lax pair given by 

i>t n =B n (><p), ip tm =B m (><p). (6) 

The t n evolutions of eigenfunction ip and adjoint eigenfunction <fr read 

rh n =B n W, <j> tn =-B* n {4>). (7) 

For n = 2, m = 1, ([5|) gives rise to the DKP equation^ 

A(/ 0t2 + 2/ 0tl - 2/o/otJ = (A + 2)/ 0tltl . (8) 

It is known that the squared eigenfunction symmetry constraint given bv[10j 

N 

B k =B k + ^2^ i A- 1 ( j> i 

i=l 

VVn = B n (lpi), 4>i,t n = - B n(<}>i), i = 1) ' ' ' 7 AT, 

is compatible with DKP hierarchy. Here iV is an arbitrary natural number, ipi 
and 0i are N different eigenfunctions and adjoint cigenfunctions of the equations 
(9c). This compatibility enables us to construct a new extended discrete KP 
hierarchy (exDKPH) as 

L u = [B n ,L], (9a) 

N 

L Tk = [B fc +X>iA-Vi,Z'], (9b) 
»=i 

fatn = B n (if>i), (t>i, tn = -B'fa), i = !,■■■ ,N. (9c) 
We have the following lemma. 
Lemma 1. Let Q = aA k , k > 1, then 

(A"VQ)- = A- x Q*(0) (10a) 
[B n ,M~V]- =B n (i))A- 1 c/>-il>A- 1 B*((t>). (10b) 
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Proof. Using /A = AT^if) - A(T- l (f)), A* = -AT' 1 , we have 

(A"VaA fc )_ = (A-^r-^^A^ 1 - A- 1 A(r- 1 (0a))A fe ^ 1 )_ 

= -(A- 1 A(r- 1 (0a))A A; - 1 )_ 
= (-l)*A- 1 A fc (T- fc (0a)) = A^A^a) = A^Q*^) 
which yields to (fTUa|) and (|10b|) . □ 

Proposition 1. The commutativity of ^9a\) and l(9b\) under \9c\) gives rise to 
the following zero- curvature representation for exDKPH 

N N 

B n ,r k - (B k +^ViA-Vi)t„ + [B n ,B k + ^ViA"Vi] = 0, (Ha) 

i=l i=l 

^i, tn = S n (^i), &,t„ = -B* n {<j>i), i = 1, 2, • ■ ■ ,N, (lib) 
with the Lax representation given by 

N 

%„ = B n (V), * Tfc = (B k + J] ^A~ VOW- (12) 

i=i 

Proof. For convenience, we omit J^. By ([5]) and Lemma 1, we have 
B n , rk = {L n Tk )+ = [B fc +VA-V,^"]+ = [B fc +^A-V,i+] + + [S fc +V'A- 1 0,i"] + 
= [S fc + VA"V, L»] - + VA"V, + [Sfc. = [B* + VA"V, L+] 
-ty>A~ V, + [fl n , L k ]+ = [B k + V»A-V, Sn] + + ^A"V)*„- 

□ 

Remark. The exDKPH (11) extends the DKPH (5) by containing two time 
series {t n } and {rfc} and more components ipi and (pi, i = 1, • • • , TV. 

Example 1. The first type of DKPSCS is given by ill)) with n= 1, k = 2 

N 

A(/ 0T2 + 2f otl - 2/ / otl ) = (A + 2)/ Mltl - A 2 ( 13a ) 

i=l 

i> iM =A(^)+/oVi, = -A*(&)-/ &, i = l,2,.-.,JV. (13b) 
Jte Lax representation is 

tf tl =(A + / )(tt) (14a) 

jv 

* r2 = (A 2 + (/o + / (1) )A + A(/ ) + A (1) + /i + fl +^^A"Vi)W- 

i=l 

(14b) 
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Example 2. The second type of DKPSCS is given by ill)) with n = 2, k = 1 

N 

A(/o* 2 + 2/on - 2/o/orJ = (A + 2)/ 0riTl + ^[A 2 ((/ + / _1 - 2)^^) 

+ A(^ 2) & - ^^" 2) ) + A((r + (15a) 
lfc, ta = A 2 (^:) + (/ + / (1) )A(^) + (A(/ ) + + A + /o 2 )^i, (15b) 
&, ta = -A* 2 W> 4 ) - A*((/ + - (A(/ ) + +/i + / 2 )V^. (15c) 

lis Laa; representation is 

* t2 = (A 2 + (/o + / (1) )A + A(/ ) + f[ 1] + h + / 2 )« (16a) 

N 

* Tl = (A + f + ^A"Vi) (*)• (16b) 

3 Reductions of the exDKPH 

3.1 The £ n - reduction 

The i„-reduction is given by 

L™ = B n or L n _ = 0. (17) 

Then we have 

(L n ) tn = [B ni L n ] =0, B„,*„=0. 
So L is independent of t n and we have 

B n fofc) = L n U>i) = Ki>i, Ki4>i) = A?&. (18) 
Then we can drop t n dependency from pip and obtain 

JV 

B n ,r k = P„)| +Y,i>i^~%,B n ], (19a) 

i=l 

B n y>i) = XfiJ it B* n {4>i) = A?&, i = 1, 2, • • • ,7V, (19b) 
with the Lax pair given by 

N 
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(fH)|) can be regarded as discrete (l+l)-dimensional integrable hierarchy with 
self-consistent sources. When n = 2, k = 1, ()19[) gives rise to 



AT 



2A(/ 0ri - /o/orj = (A + 2)/ 0TlTl + ^[A 2 (/ + /< 1} - 2)V 1 ^r 1 

i=l 

+ A(4 2 Ui - i>i4~ 2) ) + A ( r + (20a) 

A 2 (^) + (/o + /^ 1} )A(^) + (A(/ ) + +/i + / 2 )V^ = A 2 ^, (20b) 

A* 2 (^) + A*((/ + / (1) )^) + (A(/ ) + fP +h + / 2 )V, = A. 2 ^, (20c) 

which can be transformed to the first type of Veselov-Shabat equation[T^ with 
self-consistent sources (VSESCS). 

3.2 The Tfe- reduction 
The rfc -reduction is given bv[T0] 



N 



i=l 

By dropping dependency from (jlip , we obtain 

AT Ar TV 

(Sfc + ^^A-Vi)^ = [(Bfc + ^V'iA-V^I.-Bfc+^ViA- 1 ^], (21a) 

2 = 1 2=1 2=1 

AT 

V>i,t„ = (B k +Y, ^iA-Vi)l(V'i), (21b) 

i=l 
AT 

&,t„ ^-(Bfc + ^^A- 1 ^)!*^), » = 1,2,... J JV > (21c) 
i=l 

which is the k-constrained DKP hierarchy. When n = 1, fc = 2, (I21|) leads to 

AT 

2A(/ otl - / / otl ) = (A + 2)/ otltl + A 2 (22a) 



=A(VO+/olfc, 0i,ti = -A*(0i) - Mi, i = l,2,---,N, (22b) 
which can be transformed to the second type of VSESCS. 

4 Dressing approach for exDKPH 

4.1 Dressing approach for discrete KP hierarchy 

We first briefly recall the dressing approach for DKPH [11] . Assume that oper- 
ator L of DKPH HJ) can be written as a dressing form 

L = WAW~ 1 , (23) 



G 



W = A N + w^- 1 + w 2 A N - 2 + ---+w N . 
It is known [T^] that if W satisfies 

W tn = -L n _W : (24) 

then L satisfies It is easy to check the following Lemma. 

Lemma 2. If h tn = A n (h), W satisfies then ip = W(h) satisfies i.e. 

= B n (i/>). (25) 

If there are N independent functions h\,...,h^ solving W(h) = i.e. 
W(hi) = 0, then wi, . . . , wn are completely determined from these hi, by solving 
the linear equation: 



(hi A(fci) 
h 2 A(h 2 ) 



A N-1 
A N-1 

\N-ll 



(hi) 



\h N A(h N ) ■ ■ ■ 
Then the operator W can be written as 





( w N \ 

WN-\ 




(A N ( hl )\ 
A N (h 2 ) 


) 






{A N (h N )J 



w 



1 



Wrd{h\, ■ ■ ■ , hisf) 



hi 


h 2 


h N 1 


A(hi) 


A(h 2 ) ■■ 


• A(h N ) A 


A w (fci) 


A N \h 2 ) ■■ 


■ A N {h N ) A N 


hi 


hi 


h N 


A(hi) 


A(h 2 ) 


A{h N ) 



(26) 



where Wrd(hi, ■■ ■ , /ijv) 

A N - 1 {hi) A N -\h 2 ) 
Proposition 2. Assume that hi satisfies 

h iitn = A n (hi), i=l,---,N 



.N-l 



(h N ) 



(27) 



W and L are constructed by 
respectively. 



and 



then W and L satisfy \21$ and O), 



Proof. Taking partial derivative dt n to the equation W(hi) — 0: 
W tn (hi) + WA n (hi) = {W tn +L n + W + L n _W){hi) 

= {W tn + L n _W)(hi) = 0, * = 1, • ■ ■ ,N, 

since L n _W = L n W — L\W = WA n — L\W , L n _W is a non-negative difference 
operator of order < N, Wt n + L 7 ^W is also of order < N. Then according to 
the difference equation's theory, Wt n + L"^W is a zero operator. □ 
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4.2 Dressing approach for exDKPH 

We now generalized the dressing approach to exDKPH @. We have 
Lemma 3. Under fiZty , ifW satisfies and 

N 



W Tk = -L k LW + tpiA^faW (28) 



then L satisfies &9a\) and \9b\) . 

Proof. It is known that L satisfies (j9"aj) . We have 

L n =W n AW~ 1 - WAW^W^W' 1 



N 



={-L k _ + ^iA~ l fa)L + L(L k L - J2 ^A- l (j>i) = [B k + Y^ TpiA^fa, L]. 

i i i—1 

□ 

This dressing operator W is constructed as follows: Let gi, g~i satisfy 

g itU = A n ( 9i ), g ltTk = A k { 9l ) (29a) 
g i>tn = A n ( 9i ), - 9l , Tk = A k ( 9i ), i= 1,...,N. (29b) 

And let hi be the linear combination of gi and g~i 

hi = gi + ai{T k )gi i = l,...,N, (30) 

with the coefficient on being a differentiable function of r^. Suppose hi, . . . , Hn 
are still linearly independent. 
Define 

/ • n;/ji j. i nff-i Wrd(r/n, ■ ■ ■ ,r)n, ■ ■ ■ ,r/ijv) . . . 

^ = -^0,0, 0i = (-i) — wrd(r fel ,...,rM >* = ^---> N ( 31 ) 

where the hat " means rule out this term from the discrete Wronskian determi- 
nant, on — ^? ± . We have 

Proposition 3. Let W be defined by and fifty, L = WAW' 1 , tpi and fa 
be given by fifl)) . then W, L, ipi, fa satisfy fif$, fifty and exDKPH (0). 

To prove it, we need several lemmas under the above assumptions. The first 
one is : 

Lemma 4. (The discrete version of Oevel and Strampp's lemma 1181) 

N 



Proof. Note that <fii, . . . , 4>m defined in ([3"Tj) satisfy the linear equation 



TV 



J2^ j Fh i )-<t> i = 8 jtN - 1 , j = 0,l,.--,JV-l (32) 

i=l 

where <5j 1 is the Kronecker's delta symbol. Using properties /A -1 = 
E^oA^^A^r/), we have 

JV N oo oo N 

hiA-^i = J2J1 A ^' _1 Aj ( r (^)) ■ & = E A ~' _1 E A "( r (^)) • & 

2=1 2=1 j=0 j=0 i = l 

AT-1 oo N 

= Art-^N-i + £ A-^- 1 ^ A'(r(fc)) ■ ^ = A^ w + o(A- N - 1 ), 

j=0 j=N i=l 

So we have 

W^/iiA- x & = 1 + (W^/iiA -1 ^.)- = 1 + W{hi)A- l ^i = 1. (33) 

2 2 2 

This complete the proof. □ 
Lemma 5. W*(<k) = 0, for i = 1, . . . , N. 
Proof. Lemma [T] implies that 

(A- 1 <f> i W)- = A- 1 W*(<t> i ). (34) 
Using Lemma [4] and (|10a|) , we have 

N N 

= (A J W~ 1 W)-i = (A^hiA-^iW)- = Q2A'{hi)A-%W)- 

2=1 2=1 

= ^A^'(^)A- 1 ^*(0 l ), ./ 0.---..V 1. 

2=1 

Solving the equations with respect to A~ 1 W*(4>i), we find A~ 1 W*(4>i) = 0. 
This implies W*(<fc) = 0. □ 

Lemma 6. T/ie operator A (f)iW is a non-negative difference operator and 

(A^&WXfy) = 1 < *.i < N - (35) 

Proof. Lemma [5] and (f34|) implies that A _1 </>jVF is a non-negative difference 
operator. We define functions Cjj = (A -1 0iW)(/ij), then A(cy) = 0jW(/ij) = 
0, which means cy does not depend on the discrete variable n. From Lemma |H 
we find that 

N 

A k (hi)dj = A k {Y,{h i A- 1 <j> i W){h j )) = A^W^WXh,) = A fc (^-), 

2=1 2 

so cy = 8ij. □ 



Proof of Proposition^ The proof of is analogous to the proof of in the 
previous section. For ([28]) . taking d Tk to the identity W(hi) = 0, using ([29]) , 
flBTi]) . the definition (J2J) and Lemma we find 



JV 



=(W Tk )(hi) + (WA k )(h l ) + ctiWfa) = (W Tk )(h t ) + (L k W)(h t ) Y, '-'A, 

3=1 

N 

=(W Tk + L k _W -^jA-^jW^hi). 
i=i 

Since the non-negative difference operator acting on hi in the last expression has 
degree < N, it can not annihilate TV independent functions unless the operator 
itself vanishes. Hence ([2"5]) is proved. Then Lemma [3] leads to (|9b[) . The first 
equation in (|9c[) is easy to be verified by a direct calculation, so it remains to 
prove the second equation in (JScJ). Firstly, we see that 



(W~ 1 ) tn = -W^WtJV' 1 = W~ x (L n - B n ) = A n W _1 - W^Bn. 
Then we substitute W^ 1 ~ E hiA^ 1 ^ to this equality at both ends, we have 

= {A n W~ l - W~ 1 B n )- = A n (h l )A- 1 & - hiA^BUfa) 
Then ^^A" 1 ^,^ = - E hA' 1 B*((j>i) implies that (5c]) holds. 

5 N-soliton solutions for exDKPH 

Using Proposition [3] we can find solutions to every equations in the exDKPH 
(0). Let us illustrate it by solving (13]) and (15]). For (13]), let 5 t = e x ' - 1, 
n.j = e Mi — 1, we take the solution of (l29l) as follows 



g l := cxp(l\i + Siti + S^t 2 ) = , gi := cxp(/^ l + Kih + k^t 2 ) = e v ' 

hi := gi + ai(T 2 )g~i = 2^/75" exp( ^ ^ — ) cosh(f2;), fit = i(& - rft - lna ?; ). (36) 

Since X = \VAW~ 1 = A + f + fiA' 1 + • • • , we have 

f = ResAiWAW^A- 1 ) (37) 

where W is given by (2"6"]) and ([56]) . then / , Vi an d </>i given by (|3"T]) gives rise 
to the N-soliton solution for (T3"]) . 

For example, we obtain 1-soliton solution for p3|) with TV = 1 as follows 



_^ ,Ai + hi /cosh(f^i + 2#i) cosh(fii + ^i)^ Q _ Ai — fix 
~ CXP ^ 2 ^ V cosh(fii +6>i) coshfii 
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d /aT , fi+7Ji -(A!+ Ml )/2 

The 2-soliton solution ol (|T3]) with iV = 2 is given by 
/o = -A( U!l ) = ( e Al +e A2 )A(^), 



with 



r( 1 + Q2eX2 e-eO h = r 1 + aieXl e-fe 

( e A x _ e A 2 ) w y' \^( e A 2 _ e A 1 ) t; 



= /i2 pAl 



e' 



X1+X2 



w = 1 + ai— r r— e Xl + Q!2— ; T-e X2 + ona 2 — ; r— e 

g A i e 2 e 6 2 c 

e 2A 2 _ e 2 Ml e 2 M2 _ e 2Ai e 2 M2 _ g2^i 
Wl = 1 + «i ; ; — e xi + tt2 ; ; — e X2 + ai«2 1 5 — e X1+X2 . 

e 2 — e 1 6 2 — 6 1 6 2 — 6 1 

It can be shown that the interaction between the two solutions is elastic. 
For (fT5|) . we take the solution of (|29|) as follows 

gi := exp(l\i + S^n + <5fi 2 ) = e ?1 , ft := exp(Z/ij + Kin + «?t 2 ) = e m 
hi := ft + ai(Ti)g, = 2|/Oi exp( — — ) cosh(fij). 

Then 

/o = ^esA^AW^A" 1 ), /i = iiesACW^A^" 1 ) 
together with ipi and (/>i given by (|3 1 [) presents the N-soliton solution for (fT5|) . 
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